Transient dynamics of a molecular quantum dot with a vibrational degree of freedom 
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We investigate the transient effects occurring in a molecular quantum dot described by an 
Anderson-Holstein Hamiltonian which is instantly coupled to two fermionic leads biased by a finite 
voltage. In the limit of weak electron-phonon interaction, we use perturbation theory to determine 
the time-dependence of the dot population and the average current. The limit of strong coupling 
is accessed by means of a self-consistent time-dependent mean- field approximation. These comple- 
mentary approaches allow us to investigate the dynamics of the inelastic effects occurring when the 
applied bias voltage exceeds the phonon frequency and the emergence of bistability. 

PACS numbers: 73.23.-b, 72.10.Di 



I. INTRODUCTION 

In recent years, interest in the interplay of elec- 
tronic and vibrational degrees of freedom of individ- 
ual electrically contacted molecules has given rise to 
the prospering field of molecular electronicsjii^i^ii^ Single- 
molecule transistors^ and memory cells^ have been pro- 
posed and the investigation of phonon spectra of single 
molecules by an electric measurement has become possi- 
ble. The two most prominent experimental setups for 
contacting individual molecules are STM tips^i^i^ and 
mechanically controllable break junctionsiSiii. These 
were used successfully to measure the influence of vibra- 
tional degrees of freedom on current and noise properties 
of systems as distinct as hydrogen j^^i^^'^^ water— and 
benzenei^ molecules, fullereneS j^^'^^i^^ as well as carbon 
nanotubesj ^^i^^'^" 

It was found that the vibrations of the molecule can 
indeed have discernible effects on conductance proper- 
ties. One of the most striking of these is an abrupt in- 
crease or decrease (depending on the transmission) in 
the differential conductance once the applied voltage ex- 
ceeds a vibration frequency and thus allows for the ex- 
citation of a phonon. ^^'^■^'■^■^ Recently, it was shown that 
similar features at this threshold can appear in the shot 



noise properties 



This is the reason why mea- 



surements of transport properties are becoming an in- 
valuable tool for the determination of phonon spectra of 
single molecules. Such effects have already been observed 
in a number of experiments and do not require a strong 
electron-phonon coupling. 

In the opposite limit of strong electron-phonon cou- 
pling, even more spectacular effects were predicted. 
These include the suppression of sequential tunneling due 
to the Franck-Condon effecli^ or the emergence of bista- 
bility. The latter has been argued to lead to switching 
between two stable configurations of the molecule and to 
hysteresis effects. ■^'22133^^35^ 

Even from the theoretical point of view, the knowledge 
of these systems is still far from complete. Generally, the 
electron-phonon coupling on the molecule leads to the 



formation of a local polaroni^ While the Hamiltonian 
of the isolated molecule can be diagonalized in terms of 
these polarons, no exact solution is known in the presence 
of tunneling to the contacts, such that various approxi- 
mate numerical and analytical schemes have been devel- 
oped. In the case of small electron-phonon coupling, per- 
turbation theory has proved successfuli^i^'22ii2^i^ In the 
case of small tunneling, results have been obtained using 
rate equations^^Si^ perturbation theory^^ or a trunca- 
tion scheme within a Green's function (GF) formalism, 
leading to a mean-field approximationi^2i2^ 

In order to achieve a broader understanding of this sys- 
tem, and in particular of the phonon excitation process 
and the emergence of bistability, we shall investigate the 
time evolution of a molecular quantum dot with a single 
vibrational mode as a reaction to a sudden switching-on 
of the coupling to the leads. Time-dependent transport 
properties of mesoscopic systems have been examined 
in various contexts^i^ and, in particular, such switch- 
ing effects have been investigated previously for the case 
of a noninteracting quantum dot,**^'^^ and for a weakly 
Coulomb interacting quantum dot.^^ Strongly interact- 
ing systems in the Kondo regime have been investigated 
in Refs. fiTllislligllsoj . Generally, it was shown that the 
time evolution is governed by the tunneling amplitude 
but non-adiabatic effects depending on various other sys- 
tem parameters were also predicted. Experimental mea- 
surements of such transients have been achieved using 
quantum point contact setupsi^ 

The structure of this paper is as follows: In Section 
II, we shall introduce the model used to describe the 
molecular quantum dot, the coupling to the phonon and 
the implementation of the switching. In Section III, the 
results for the case of weak electron-phonon interaction 
will be derived using perturbation theory. In Section IV, 
we go to the opposite limit of strong electron-phonon 
interaction and present results obtained by means of 
a self-consistent time-dependent mean-field approxima- 
tion. The two results will be compared in the appropri- 
ate limits and our conclusions will be detailed in Section 
V. 
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II. SYSTEM 

In general, the typical systems investigated in the field 
of molecular electronics can be rather complicated: the 
Coulomb interaction between electrons on the dot may 
be quite strong, complicated molecules may support a 
whole spectrum of phonon excitations and the transport 
of electrons to the leads may occur via multiple channels 
of different transparencies. 

In order to avoid these complications, we focus on the 
archetype model which is able to capture the most promi- 
nent effects in the physics of molecular junctions. It is 
described by the Anderson-Holstein Hamiltonian 



H 



E 



+ H, 



dot 



ph 



cl-ph 



Hi 



(1) 



=L.,R 



and it has been widely used in the literaturcj^ i^^'^'^ The 
left and right lead Hamiltonians Hlji describe non- 
interacting electron gases with the chemical potentials 
/ia, which can be tuned by the applied bias voltage 
V — — ^J'R, (we use units where e = h = 1) 



k 



(2) 



Without loss of generality, we shall assume symmetric 
bias voltage ^il,r = ±V/2 in the following. The sec- 
ond term in the Hamiltonian ([T]) represents a single level 
quantum dot with energy A, 



H, 



dot 



(3) 



which can correspond to either the LUMO or HOMO en- 
ergy level of the molecule depending on the ground state 
population. For simplicity, we have discarded the elec- 
tron spin. The interaction of the electron on the dot 
with a Holstein phonon of frequency fi, described by 
Hph — fia^a, is included up to linear order in the os- 
cillator displacement q ^ a-\- a) via the electron-phonon 
coupling parameter A 



el-ph 



= \{a + a^)dU . 



(4) 



The connection between the dot and the leads is de- 
scribed by a local tunneling Hamiltonian with amplitudes 
1L,R- In order to observe switching effects, we allow for 
time-dependent amplitudes and thus use 



HT{t) = ^-1^{t)Ui{x^Q)d + \,.c 



(5) 



where the electron field operator ipaix) is the Fourier 
transform of the electron annihilation operator Ck,a- The 
assumption of local tunneling at a: = allows us to 
discard the spatial variable x, ie. we shall use tpL,R = 
i^L,R{x = 0). As to the time dependence, we shall 
assume that the tunneling is instantly switched on at 
i = which means ^a{t) = 7a&{t), where 9{t) denotes 



the Heaviside step-function. In order to simplify the cal- 
culation, we shall assume a spatially symmetric system, 
i.e. 7L — JR = 'J and restrict ourselves to zero tempera- 
ture in the following. 

We are mainly interested in the time evolution of the 
dot population n{t) = d^{t)d{t) and the time-dependent 
current through the left/right contact. Using the Heisen- 
berg equation of motion, the current operator {a = 
L, R) can be calculated as the time derivative of the total 
charge Qa in the respective lead. 



Il,r = T 



dQL,R 

dt 



±^l{ARd"dHL,R) ■ (6) 



The signs are chosen such that current flowing from left 
to right will be positive. The expectation values Ia{t) = 
{I ait)) and n{t) = {n{t)) can be expressed in terms of 
the dot Keldysh GF, 



D{t,t')^-t{Tcdit)dHt')) 



(7) 



where Tc denotes the time-ordering operator on the 
Keldysh contour C and the time variables t and t' can be 
situated on the forward or backward part of the Keldysh 
contour C± , giving rise to the GF matrix. 



D 



D- 



D- 



(8) 



For our purposes it turns out to be more convenient to 
introduce the retarded and advanced GFs, which are de- 
fined by 



D- 



D-+ 

D++ 



(9) 



Now, the time-dependent dot occupation can be ex- 
pressed as 



i{t) = -iD~+{t,t) , 



(10) 



while the current can be written as a sum of two contri- 
butions, la = I'a + I'a, which are defined as^ 

/•DO 

I^Rit) = T7' Re / dh glR(t,h) D^'itut) , (11) 
Jo 

I'Uit) = ±7^ Re / dt, Z?«(i,ti) glRih.t) . (12) 
Jo 

Here, ga{t,t') — —i{Tc ipa{t)ipl^{t'))o denotes the uncou- 
pled lead GF, where the expectation value is taken with 
respect to the ground state of the uncoupled Hamiltonian 
Hl + Hji. It is given by 



laii^) = i2Trp{u) 



/a - 1/2 U 
-(l-/o) /a -1/2 



(13) 



which depends on the density of states p{i-u) in the leads 
and the Fermi functions fa{^) — npiuj — p,a). The super- 
script "/C" in Eqs. (dJ) and (HH) denotes Keldysh GFs 
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which are defined by = D-++D+- = 20'+ + D^- 
and analogously for . The currents through the 
two contacts are not independent since the symmetric 
bias entails Il{V) = -Ir{-V). 

In the noninteracting case (A = 0), the GF DQ{t, t') can 
be calculated analytically even for the time-dependent 
tunneling in Eq. ([5]) . The calculation is greatly simplified 
by assuming a constant density of states po in both leads. 
Introducing the contact tunnelling rate F = 277/907", o^^6 
findsi^iS for i,t' > 0, 

j^im^t - t') = -ie{t ~~ t')9{t') e-*^(*-*') e-r(t-t') 
D'^°^^{t - t') = t9{t' - t)0{t) e-»^(*-*') e^(*-*') (14) 



and 



a=L.R' 



Uuj + A) 

F2 +Cj2 



X (e— *-e-")(e'-*'-e-"') . (15) 

This leads to the following evolution of the dot occu- 
pation number in the case of an initially empty dot 
n(0) = 0, 



(t) = l6l(t)e-r*V / diu 
X [cosh(Fi) — cos{ujt)] . 



(16) 



The time-dependent current can be split into a displace- 
ment current /disp(0 = ^l(0 ^ Init), which reflects the 
change in dot population, /disp(i) = —nit), and the to- 
tal current I{t) = [lL{t) + lR{t)]/2 which measures the 
charge transported through the system. Both current 
terms as well as the dot population ri*^°)(t) where cal- 
culated in [4^ for the non-interacting case. Therefore, 
in the following sections, we will focus on how these ex- 
pressions change due to the presence of electron-phonon 
interaction. 



III. PERTURBATION THEORY IN 
ELECTRON-PHONON COUPLING 

In the limit of weak electron-phonon interaction, we 
treat the system perturbatively in A. The two non- 



vanishing corrections up to order A^ are commonly re- 
ferred to as the "tadpole" (subscript 1) and the "rain- 
bow" term (subscript 2). The former represents the in- 
teraction of a phonon with the electron density, whereas 
the latter describes virtual phonon mode exitations. The 
corresponding dot GF contributions are (superscripts (2) 
denote expressions to second order in A) 

i?f^(s,s') = 2A2 f dsi f ds2 D^"\s,si) (17) 
Jc Jc 



D'i\s,s') = 2i\^ / dsi / ds2 D^°\s,si) (18) 



and 



Xl?(")(si,52)F(5i-S2)i?(°)(s2,s') 

The time integrations run along the Keldysh contour C 
(where all time variables are on either the forward or 
backward part of the time loop) and F{t — t') denotes 
the unperturbed Keldysh GF of the phonon. In the case 
of zero temperature, it is given by 



g-iO|t-t'| gin(t-t') 
g-iO(t-t') gif2|t-t'| 



(19) 



Both diagrams translate into time-dependent contribu- 
tions to the dot occupation number which are given by 
n^i\{t) = —iD^^l ^{t,t). Up to the second order in the 
electron-phonon interaction, the total dot population is 
given by n{t) — n*^"-* -|- n^'' -|- n'^^ One can easily derive 
the following expressions for the second order contribu- 
tions. For the tadpole term, one finds 



np^(t) = 4A2 Im 



dh i?(")^(t-ti)i?(°)-+(ti,i) 

(20) 



X / dt2 F''{h~t2)n(^>it2) 
Jo 



where we used the retarded phonon GF 
F^{t - f) -ie{t - t') s\-a{n{t ~t')]. Moreover, the 
contribution from the rainbow term becomes 



^('Ht) = 4A2Re / dh / dt2i?("^''(i-ii) i?^°^"+(ti, i2)F^(ti - ^2) + i?^°^''(ii - t2)^^""(ii - ^2) D^'>'^-+{t2,t) 



(21) 



"'0 

t 



+ 2X^ dh dt2D^"^''it,h)D^°')~+{h,t2)F-+{h-t2)D^"'>^{t2,t) 
Jo Jo 



r 



In order to calculate the time-dependent current, we the second order current as la — la + la ■ Using 
make use of the general Eqs. (fTTjl and (fT2| and write 
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FIG. 1; (Color online) Time-dependent current I^^\t) for fast 
phonon mode f2 ^ F. A retardation of order between 
currents across left and right contact leads to the emergence 
of plateaus in the time trace of the total current. Parameters 
are V = O, A = n/40, A = 0. 



FIG. 2: (Color online) Time-dependent current I^^\t) for 
slow phonon mode f2 ^ F. Shortly after switching on the 
tunneling, the current overshoots to values beyond the steady- 
state value. The overshooting amplitude is characterized by 
AI{V). Parameters are n = O.IF, A = n/40, A = 0. 



iTTpo6{t), one 



the unperturbed retarded GF nit) 
finds for the first term I'j^ ji{t) = ±9 j ^ 
the corresponding second order contribution reads 



±dit)^[l - 2n{t)], so 



(2) 



it) 



,(2) 



it)] 



(22) 



Due to its symmetry, this term does not create any net 
current I{t) but only contributes to the displacement cur- 
rent /disp(i)- The second contribution from Eq. is 
given by 



Ji'(^)(i) = Tr Im 




(2)Ri 



27r 



-*)[2/i,fl,(c.)-l] 



(23) 



For the evaluation of these observables, we shall focus 
on the particle-hole symmetric case, A = 0. Deviations 
from this point were investigated in detail in Ref. [46l| and 
it was shown that they can have a clear infiuence on the 
time-dependent current and dot population as they cause 
oscillatory behavior on a time scale A~^. As these oscil- 
lations might obscure features induced by the coupling 
to the phonon, we choose to set A = 0. The remaining 
integrals can easily be solved by a simple numerical in- 
tegration which eventually leads to the time-dependent 
current and dot occupation. 

For a fast phonon mode, in the limit F <C i^, the to- 
tal time-dependent current I^"^' = [/{^ + /jj ]/2 exhibits 
plateaus (as depicted in Fig. [1]) that originate from a 
retardation between left and right current, /jf''(i) and 
/^^^ (t) . Since both the width of the plateaus and, natu- 
rally, their period are proportional to f^"^, we attribute 
this to an effect comparable to electron shuttlingj^Si^i^ 
in this particular case due to the electron-phonon inter- 
action. 



When going to the opposite regime of low phonon fre- 
quency, F ^ ri, there is a significant overshoot of the 
total current correction for short times (see Fig. [2]). Scan- 
ning through different voltages, one sees that the negative 
peak of the perturbative current grows linearly with V. 
We can characterise the relative strength of the overshoot 
by the difference between the (negative) current peak and 
its stationary value as a function of voltage 



AI{V) = Hm [/ 



(2) 



(<)] - min[/ 



(2) 



(24) 



which results in Fig. [3l The relative overshoot still in- 
creases linearly ioi V < fl but begins to decrease there- 
after. We know from the stationary-state calculation^^ 
that a,t V — Q, the phonon mode can be excited and 
inelastic processes set in which in our parameter range 
(large transmission) lead to a decrease in the station- 
ary current. This suggests the conclusion that the time- 
dependent current can also be split into an elastic part 
(which is continuous at ^ = 51) and an inelastic one 
(which vanishes identically for V < ^) and that the over- 
shooting is dominated by the elastic part. 

We would like to point out that in the first case of the 
fast phonon regime, the emergent plateaus are damped 
on the timescale of F~^. In the F ^ 51 limit on the 
other hand, the timescale for the dominant effect, the 
overshooting, is also of the order of F~^. Therefore, the 
plateau feature, which is dominant in the former regime, 
is not visible anymore, because the time evolution of the 
phonons is too slow. 



IV. MEAN-FIELD APPROXIMATION 

The perturbation theory is valid only for low electron- 
phonon coupling A. In order to describe features like 
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n/r = 0.1 

Q/r = 0.05 
Q/r = 0.075 




FIG. 3: (Color online) Overshooting strength AI as a function 
of bias voltage V. The overshooting displays a maximum at 
the threshold for inelastic processes. 



where the electron density is denoted as n = (d^d). 
Hence, we see that the electron-phonon coupling leads 
to a shift of the rest position to 



{q) = - 



V2Xr. 
O 



(28) 



In the adiabatic case, this equihbrium position follows 
the time evolution of n{t) and thus changes on the time 
scales and r~^. If the electronic dynamics is slower 
than the phonon dynamics, the physical picture is that of 
a very fast oscillator that slowly adjusts its equilibrium 
position. Inserting this shift into Eq. (P5|) . wc find that 
the back-action from the mean displacement leads to a 
time-dependent dot energy 



o\2 

A'(t) = A-— n(0. 



(29) 



bistability, the electron-phonon coupling strength needs 
to be much higher. A possible route to a solution for 
strong interaction is provided by a Hartree-like mean- 
field ansatz, which we would like to justify from the fol- 
lowing physical point of view. 

We know from the non-interacting case in Eq. (jl6p that 
there are basically two time scales governing the time- 
evolution for instant switching-on of the tunneling. It is 
the exponential saturation on the scale as well as an 
oscillation with a frequency of order A. If the oscillator 
period Sl^^ is much shorter than any of these time scales 
(il 3> A,r) then the change of the dot population will 
look very slow from the oscillator's point of view, and the 
oscillator rest position can adapt adiabatically to the mo- 
mentary n{t). Therefore, the mean-field approximation 
leads to the effective Hamiltonian for the electrons, 

Hei= Ha + Ht + AdU + X{a + a^) dU (25) 

a=L,R 

and the phonon Hamiltonian can be written as 

Hph = fla''a + \{dU){a + a''). (26) 

In terms of the dimensionless position and momentum 
operator q = -^[a + a^) and p — ^{a^ ~ a), we can 
rewrite the phonon Hamiltonian as 

Hph = § + P') + ^/2Anq - ^ (27) 



After this motivation, we shall derive the real-time dy- 
namics of the system in this regime in a more formal man- 
ner. The diagrammatic access to the mean- field ansatz 
is gained by a summation of tadpole terms, which is jus- 
tifiable for large fl, since the perturbation theory reveals 
that the tadpole term is the dominant one in this case. 
For this purpose, we start from Eq. (fT7|) and replace two 
dot GFs on the right hand side by their exact counter- 
parts. Thus, we obtain the Dyson equation, 

Dm{is,s')^D^°'>is,s') + 2\^ [ dsi [ ds2 D^°\s,si) 

Jc Jc 

X Dmfisi,s')F{si - S2)»T-mf(s2) (30) 

This can be written in the conventional form D = D'^*'^ + 
D'^'^^SD by defining a self-energy in Keldysh space. For 
short phonon periods and coherence times compared to 
and A~^, it can be approximated as 



Smf(si,S2) = — ^rimf(si)<5(si - S2) (31) 

This time-dependent self-energy corresponds to a time- 
dependent dot energy level as in Eq. We can solve 
this mean-field Dyson equation and end up with the fol- 
lowing self-consistency equation for the dot population 



duj 



E„/a(^-A)-2no / A 



Y2 



(32) 



where n^^^ {t) is given in Eq. (flBl) . no G {0, 1} is the initial dot population and 



-N{t) 



h{t,uj) = / ds g{s,uj)nn,f{s) 
Jo 



(33) 
(34) 
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,(2) 



it) 



,(2) 



(t) very well. However, no bistability ex- 



FIG. 4: (Color online) Time-dependent dot population n{t) 
for initially empty (solid lines) and occupied (dashed lines) 
dot level. In the bistable regime, the different initial states 
evolve into different stationary states. In the inset: Bistable 
regions as a function of dot detuning A. The parameters used 
are r = 1, A = 6 and f2 = 10. 



and N{t) = ds nnif(s). Note that Eq. ((5^ incorpo- 
rates the memory behavior of the mean-field ansatz since 
the electron density at time t directly depends on all pre- 
vious values after the tunnel switching at t = 0. This 
retardation makes this problem solvable numerically by 
a discretisation of the time axis. 

The resulting dot population in the mean-field approx- 
imation is shown in Fig. U) The locations of bistability 
as a function of A agree with those found in the self- 
consistent Hartree approximation of Ref. [s^ . Both sta- 
ble solutions can be reached depending on the initial dot 
population. For uq — 0, the lower state is reached, while 
for no = 1, the system evolves into the stable state with 
the higher dot population. 

In the limit of weak electron-phonon coupling, we com- 
pared the mean-field results to the perturbative results 
of the previous section. It turns out that for small A, 
rim fit) coincides with the function resulting from the 
tadpole term, n'"^\t), given in Eq. (PO)) . This is not sur- 
prising since the mean-field approach corresponds to a 
summation of all diagram of the tadpole type. Since for 
a fast phonon, the tadpole term is the dominant one, 
the mean-field result reproduces the perturbative result 



ists in this perturbative regime. Conversely, perturbation 
theory is not applicable in the bistable regime. 

Qualitatively, the function n(t) as calculated by the 
mean-field approach retains the central characteris- 
tics known from the non-interacting case. Its time- 
dependence is governed by an exponential growth on 
a time scale with superimposed oscillations on a 
time scale A^^. However, the steady-state value can be 
strongly influenced by the presence of the phonon mode. 
V. CONCLUSION 

We investigated the transient effects occurring in a 
molecular quantum dot coupled to a single phonon mode 
when the tunnelling to the leads is switched on instantly. 
In the case of small electron-phonon coupling, we used 
perturbation theory to calculate the time-dependent dot 
occupation n{t) as well as the time-dependent currents 
through the two contacts /L,i?(i)- 

The sudden switching leads to a number of non- 
adiabatic effects. We found that in the regime il 3> F, 
a retardation between /^(i) and I nit) leads to the emer- 
gence of plateau structures in the time trace of the total 
current, which we attribute to a feature similar to elec- 
tron shuttling. 

In the opposite limit of a slow phonon, il ^ F, we find 
a voltage-dependent overshooting of the current com- 
pared to its steady-state value. The overshooting has a 
maximum at the threshold voltage for inelastic processes. 

In the case of stronger electron-phonon coupling, we 
used a generalisation of the mean-field approximation to 
investigate n{t). We found that this scheme correctly 
reproduces the known static results and the bistabil- 
ity emerges naturally as the two initial dot occupations 
evolve into different stationary states for certain param- 
eter constellations. 
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